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Abstract. Finite difference approximations to multi-asset American put option price 
are considered. The assets are modelled as a multi-dimensional diffusion process with 
variable drift and volatility. Approximation error of order one quarter with respect to the 
time discretisation parameter and one half with respect to the space discretisation parame- 
ter is proved by reformulating the corresponding optimal stopping problem as a solution of 
a degenerate Hamilton- Jacobi-Bellman equation. Furthermore, the error arising from re- 
stricting the discrete problem to a finite giid by reducing the original problem to a bounded 
domain is estimated. 



1. Introduction 

American put option is a derivative contract based on the price of some asset, denoted 
by St, which evolves with time. At time t — Q, when the contract is entered, a reference 
asset is chosen together with "strike" i^strike and expiry time T. The contract gives the 
holder the right (but not the obligation) to sell the asset at any time T* e [0, T] for the 
amount i^smke- Thus, when exercised, the payoff from the American put is K^mke — S't* ■ 
The option will only ever be exercised by a "rational" investor if the payoff is positive, so 
the payoff from American put is [if strike ^ St']+, where [x]^ denotes the positive part of 
any real number x. Even in the classical Black-Scholes model, there is no known formula 
for the price of an American put with a finite exercise time (there are formulae for prices 
of infinite exercise time American put, American call option and European put and call 
options). A variety of numerical methods and approximations for the American put option 
price have been developed over the years. An overview of the various methods can be 
found for example in Barone-Adesi [4|. 

Four main approaches for calculating the put option price can be identified: The first 
attempts to find formulae that give results close to the real price. These give fast approx- 
imations but the accuracy cannot be simply improved upon by doing more computations. 
The second approach approximates the evolution of the underlying asset with a recombin- 
ing tree (typically binomial or trinomial). Backward induction on the tree then yields the 
American put price. Error estimates for the binomial tree approach are proved in Lam- 
berton |21| and improved in Lamberton and Rogers ll24l and Lamberton ll22l . The third 
approach is based on Monte Carlo methods. Finally the fourth approach relies on refor- 
mulating the option price as a solution to a partial differential inequahty or a nonlinear 
partial differential equation. This equation can then be discretised using a variety of meth- 
ods. For the finite element method see for example Allegretto et al. |T| and also Pironneau 
and Achdou [25. Section 6.4]. Finite volume methods have been used by Angermann and 
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Wang ['2| and Berton and Eymard fS). Nevertheless finite difference methods are particu- 
larly popular and it is those that we focus on here. 

If a diffusion process is used to model the asset price, then it has been shown that the 
American put price is the payoff function of an optimal stopping problem. See for exam- 
ple Shiryaev ||26| . The payoff function for this optimal stopping problem satisfies a system 
of second order partial differential inequalities. The solution to this system can be ap- 
proximated using finite difference methods. See for example Lamberton and Lapeyre |,23l 
Chapter 5] for a concise introduction. They also present a simple algorithm for computing 
the solution to the finite difference problem. A more efficient iterative method based on 
the SOR method, called projected SOR, is given in Pironneau and Achdou [25., Section 
6.4. 1]. Finding more efficient ways of computing the solution of the finite difference prob- 
lem are of considerable interest. The reader is referred to Forsyth and Vetzal flOl and Cen 
and Le |7J where the penalty method and singularity separating, implicit finite difference 
scheme, are studied and compared with other methods. 

This paper is focused on implicit finite difference approximations to the American put 
option price in the case when there is one or more underlying assets with variable diffusion 
coefficient, drift, and discounting. We prove that the error introduced by the implicit finite 
difference approximation is, under suitable regularity assumptions, of order r^/'' + h^/'^, 
where r denotes the space discretisation parameter and h denotes the space discretisation 
parameter. The only other result in this direction is Hu et. al. lfT4l . as far as the author is 
aware. In Hu et. al. [ 14 1 rate of convergence of order t^/^ + is proved and furthermore it 
is show that this is an optimal rate of convergence. However this is done in a much simpler 
setting then this article considers. In particular, only one risky asset is considered and it is 
assumed to have constant drift and diffusion coefficients. Furthermore this article adds the 
estimates for the error arising in computing the discrete problem on a finite domain. The 
reason why the same rate of convergence as in Hu et. al. llT4ll is not obtained in this article 
is that here the diffusion coefficients are allowed to degenerate. Hence the solution can 
only be Lipschitz continuous in the space variable. The result of Hu et. al. [141 requires 
more regularity and only holds in the non-degenerate case. 

The way to obtaining the rate of convergence has been paved by recent work of Krylov 
on the rates of convergence of finite difference approximations to the Bellman equation 
(also known as Hamilton-Jacobi-Bellman equation). The first rate of convergence esti- 
mates were obtained in Krylov I17J . This has later been extended to the case of variable 
coefficients in Krylov ifTSl . The rate of convergence has been further improved, in the 
case of constant coefficients, in Krylov and Dong, [[8|. Finally the convergence rate of or- 
der T^/^ + h^/^ for Bellman equations with Lipschitz continuous coefficients is obtained 
in Krylov |fT9l . These results have been extended in Gyongy and Siska lfT3][ to cover the 
Bellman equation corresponding to optimal stopping of controlled diffusion with Lipschitz 
continuous coefficients. In Krylov and Dong [9[ the results from Krylov [fT9l[ have been 
extended to allow domains not equal to W^. In Krylov [20[ one of the key ingredients of 
the proof, the discrete gradient estimate, has been generalised to allow estimates for other 
nonlinear partial differential equations. The constant coefficients case has been studied, 
also adapting Krylov's methods, in Jakobsen [15[. Krylov [19[ considers general finite 
difference schemes which have been already introduced in Bonnans and Zidani ||6| in the 
controlled Markov chain setting, however without establishing any rates of convergence. 

The problem of restricting partial differential equations arising in finance to bounded 
domains has been studied in Barles, Daher and Romano [3[. By introducing artificial 
boundary condition of either Dirichlet or Neumann type on the boundary of the ball to 
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which they restrict the domain they are able to prove convergence to the solution of the 
equation on the whole space. This applies to a class of nonlinear partial differential equa- 
tions to which the viscosity solutions exist. However only in the case of linear equa- 
tions with constant coefficients do they get an exponential rate of convergence. This paper 
demonstrates exponential rate of convergence for a nonlinear problem (the American op- 
tion price) with variable coefficients. This is proved by first obtaining a general result on 
the distribution of exit times of a diffusion process from a ball, see Lemma l4n and second 
by applying the maximum principle for the discretised equation. 

The paper is organised as follows. In Section |2] the main result is presented together 
with the assumptions. In Section [3] the partial differential equation for the American put 
option price is obtained together with the rate of convergence estimate for approximations 
on infinite grids. In Section|4]the error arising from restricting the infinite grid to a finite 
grid is estimated. 

2. Main result 

Let (ri, J^, P, {Tt)t>o) be a probability space with a right-continuous filtration, such 
that To contains all P null sets. Let {Wt,J^t) be a c?-dimensional Wiener martingale, i.e., 
let {Wt )t>o be adapted to {Tt )t>o and for all t,s > 0, Wt+s — Wt independent of Ft- 

We will consider the standard Black-Scholes model extended to several dimensions. We 
consider d risky assets and one risk-less asset. We assume that we are given p = p{t, x), a 
non-negative real valued function of i G [0, T] and x e M'', representing the continuously 
compounded zero coupon rate and also a = a{t,x), a d x d matrix valued function of 
t e [0, T] and x e M'', representing volatilities of the risky assets and the correlations 
between the risky assets. 

Assumption 2.1. The functions p and a are Borel measurable in t. There exists a positive 
constant K such that 

\a\ = {traa'^y^^ < K and p < K on [0,T] x . (1) 

We denote the risky assets Su — {S^, . . . ), where Su = {Su^)ue[t.T] is defined on 
[t, T] by the stochastic differential equation 

d 

dSl^Slp{u,Su)du + SlY.<j'^{u,Su)dW.i, Sl^S\ t^l,...,d, (2) 

where 5* > 0. It is well known (see e.g. Krylov fid', Chapter 2, Section 5]) that the sto- 
chastic differential equation has a unique solution under Assumption l2. 1 [ together with the 
assumption that the functions S i-^ Sa{u, S) and S i-^ Sp{u, S) are Lipschitz continuous 
for all u E [0, T]. We will use the notation £4 5 to denote the expectation of the expression 
following with the understanding that the relevant stochastic process is started from point 
S at time t. We consider the optimal stopping problem 

v{t,S)^ sup Et,s (e- ^^"^''"'"'"HSt')) , (3) 
T*e-x[t,T] ^ ^ 

for a given Lipschitz continuous function g. The American put option price is given by 
V. See e.g. Shiryaev [26J. In the one dimensional case g{S) :— [Ksmke — S]+ but in the 
multidimensional case one may want to consider a general payoff g. We wish to remove the 
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linear growth present in the coefficients of (HI. Hence for m G [t, T], we define := In S*^, 

d 

j=i 

and cr(i, a;) = a{t, e^), I3{t, x) = /3{t, e^) and p{t, x) = p{t, e^). By Ito's formula we get 

d 

dxl=(3\u,xl,)du+Y,'^''{u,xl,)dWi, xl^x'^\nS\ z = u e [t,T]. 

3=1 

Let g{x) := g{e'^' ) and x = In 5. Then the option value v{t, S) given by ([D is equal to 

w{t,x)= sup Et.:,(e--''o^*''("'^")''"g(xT.))- (5) 

T'e'I[0,T-t] 

Let 77 = ?7(x) be a smooth function and define 

d d 
ij — 1 2—1 

Assumption 2.2. There exist a natural number di, vectors i'fc G R'* and functions 

afe : [0, T] X M'' ^ M, 6fe : [0, T] x R'^ R, fc = ±1, . . . , ±di, 

such that £k = —^-k, Kfe| < a*: = a-k, bk > for A; = ±1, . . . , ±di and such that 

^iaa^y^it,x) = afe(i,a;Kfe£^„ and pHt,x) = 6fe(i,x)f^ 

for all k € {±1, . . . , ±di}, i, j € {1, ... ,4 and {t,x) e [0,7] x R''. 

Let and denote the first and second derivatives in the direction of a vector £ in 
R''. Notice that under this assumption L given by ^ satisfies 

k=±l,...,±di 

Remark 2.3. While this assumption may appear restrictive, it turns out that this can be 
satisfied for any operator given by (|6]), such that: 

(1) We can find Lhri{x) = J2yeB Phiy)v{^ + ^J/) with a finite S C R'^ such that 
span B = 

(2) WehaveL77(0) = lim,,^o U7?(0). 

(3) If 77 has a strict maximum at then Lft77(0) < 0. 

For a proof see Dong and Krylov IH Section 3]. In the one dimensional case the con- 
struction of Ofc and bk is straightforward. In several dimensions in the case when ^crcr^ is 
diagonally dominant see e.g. Siska |28, Example 5.2.5]. See Bonnans and Zidani |61 for 
the case when ^crc^ is not diagonally dominant. 

We will need the following regularity assumptions. 

Assumption 2.4. For ip e {ct, y/ok, bk, pfor k = ±1, . . . , ±di} there exists K > such 
that for all t, s e [0, T) and for all x, y G R"^ 

\i}{t,x)-i^{s,x)\<K\t-s\''^, \il^{t,x)\<K. and \^{t,x)~i:{t,y)\ < K\x-y\. 
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Notice that Assumptions 12.21 and 12.41 imply Assumption 12.11 Let ft, > 0, t > and 
£ £ M.'^. We define TT{t) := min(r, T — t). So the time step is fixed except the case 
t £ (r, T) when T — tis used instead of r. Let 



rT„/^„^. il{t + TT(t),x)--q{t,x) rj{t , x + M) - rj{t , x) 



5^ r]{t,x) := ■ ■ — , 6hjri{t,x) := 

T 

^hj V '■= - Sh,i 5h-i V = j^i^hj V + 5h-i v) 



(8) 



fort e [0,T),a; e R''. Let 

fc=±l,...,±di 

where Uk and bk are the functions from Assumption l2.2l 



Remark 2.5. Contrary to the usual finite difference approach we have not yet introduced 
any grid on which the solution to the discrete problem is defined. Typically, one first 
introduces a grid and then the finite difference operators acting on functions on the grid. 
Here the opposite approach is taken. First the finite difference operators are defined for any 
point {t, x). Thus we will obtain a collection of disjoint problems, each centered around 
an arbitrary (to, xq) G [0, T] x M"* and solved on the grid Mt Mt H ([0, T) x M'^), 
where 

Mt ■■= {{t, x) e [0, T] X M'' : {t, x) = ((to + jr) A T, 

xo + Knii + ■■■ + idjd,)),j e {0} u N, ifc e Z, = ±1, . . . , ±di}. 

The discrete problem to be solved in order to approximate the price of the American put 
option is 

S'^Wr,h+l^hWrM,9-Wr,h = On Q, ^^^^ 

Wrji = 9 on Mt \ Q, 

where Q C Mt- The solution Wr^h to dTol ) can be defined for any point in [0, T) x W^, 
since the grid Mt can be centered arbitrarily. 

Remark 2.6. It is worth noting the nonlinear structure of the above partial differential 
equation. If the equation was linear and non-degenerate, i.e. if, for example, on the left 
hand side we only had the first term of the maximum and non-degenerate L, we would 
be in the standard situation of linear parabolic equations and we would immediately know 
that what the rate of convergence is, from, for example, Thomee 1271 . 

Finally, we consider the localisation error We will use R > Ri > R2 > 0. Let 
Br = {x £ M."^ : \x\ < R}. We will need to solve the discrete problem on a grid in 
Bn. The radius Ri is used for introducing the artificial boundary conditions, while in Bti^ 
we obtain the desired estimate. Let be a function that is equal to g inside zero 
outside i3j?i+i and Lipschitz continuous. Let Qji := ([0,T) x Bb) H Mt- The discrete 
problem that needs to be solved is 



: on Qr, 

(11) 

= fffli on Mt \ Qr- 

The following theorem is the main result of this paper. We will always use C > to denote 
a generic constant that is independent of T,h,Ri, R2 and R. 
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Theorem 2.7. Let Assumptions \2?2\ and \2^ be satisfied. Then the system jlOl) has a unique 
solution for any Q <Z Mt and for any g — g{x) which is bounded and Lipschitz continu- 
ous in X. Furthermore, let w be given by (|5]l. Let ^ denote the solution to (II 11 1. Then 
there are constants fi > and 7 € (0, 1) such that 

\yj _ ^^f^ I < C(e-pfl?+fl^/2 + ^1/4 ^ ^ [0, T] X i?;^,. 

The proof will be given in Section ID An outline is given below. We start by using 
randomised stopping to express the payoff function of the optimal stopping problem as a 
payoff of an optimal control problem with unbounded reward and discounting functions. 
The payoff of this optimal control problem then corresponds to the solution of a normalised 
Bellman equation which we use to derive the finite difference approximation. Adapting 
results from Gyongy and Siska IT31 gives the first main result of this paper, the rate of 
convergence of order t^/^ + h^^"^ for a grid on the whole space. We then prove an estimate 
the probability that a stochastic process exits a certain ball before time T and use this 
together with a discrete comparison principle to estimate the error arising in restricting the 
approximation to a finite grid. 

3. Normalised Bellman equation 

This section appUes known results about optimal stopping, optimal control and nor- 
malised Bellman equations to estimate the rate of convergence. The following theorem is 
a special case of Gyongy and Siska LI 2. Theorem 3.2]. It is also proved in Krylov lfT6l 
Chapter 3]. 

Theorem 3.1. Let Din contain all progressively measurable, locally integrable processes 
f — ii^t)t>a taking values in [0, n], such that rtdt = 00. Let SH = UneN 9 
a Lipschitz continuous function of X. Then for all (t,x) S [0,T] x M.'^, 

w{t,x) = supEt.^ ( / r,5(x,)e-^*''''("^^"'+''"''"ds + <7(a;T)e-^*^''("^^")+''"''"] . 

Remark 3.2. Let L be the differential operator given by (|6]l. Krylov lfT6l Theorem 6.3.3] 
proves that w is the unique solution of the following normalised Bellman equation 

sup(j^{wt+Lw) + j^ig-w))=0 on[0,r)xM^ (12) 

w{T,x)=g{x) for x eR'K (13) 

Let e — and use this in (fT2] |. The supremum is now taken over e G [0, 1] and hence 

sup (e {wt + Lw) + {l-e){g- w)) = on [0, T) x M'^, 
ee[o,i] 

Noticing that for any real numbers p and q, 

sup {ep + (1 — s)q) = max(p, q) 

we obtain that (fT2] l is equivalent to 

max [wt + Lw, g - w] = on [0, T) x R''. 

The reader will immediately recognise that this is exactly the nonlinear equation of which 
([Tol l is the finite difference approximation. 
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Lemma 3.3. Let Q be a subset of Air and g be a bounded and Lipschitz continuous 
function. Let Assumptions \2. 71 \2.2\ and \2.4\ be satisfied. Then (llOl l has a unique solution. 

Proof. We wish to transform (fTol i into the same form as it appears in Gyongy and Siska ifTsll . 
Using the same argument as in Remark 13.21 we can see that the first equation in ( fTOl l is 
equivalent to 

sup (s^ Wr.h + l^hWr.h - rwr.h + rg) =0 on Q. (14) 

re[o,oo) l + r \ ■ J 

Now we just need to check that the assumptions required by Gyongy and Siska ITSl Theo- 
rem 3.4] are satisfied in our case. The space of control parameters is [0, oo). The normalis- 
ing factor m is 1/(1 -I- r). Assumption 3.1 of flTl is satisfied due to Assumption |2.2| The 
discounting function c is p + r and the reward function / is rg. Due to Assumption |2. 1 1 we 
have 

+ \t^M + \TrAp'rr)\ + \^vg\ < C. 

Thus fTT, Assumption 3.2] is satisfied. Finally Yq-;:(1 + p + r) > 1, since p > 0. Hence 
|fT3 Assumption 3.3] is satisfied and the result follows from 1 13 , Theorem 3.4]. □ 

Theorem 3.4. Let w be the American put option price given by (O. Let Q — A4t and let 
Wr^h be the solution of dlOt . Then 

\w - wr,h\ < C(t^I^ + h^l"^) on Mt- 

Proof. As before we consider the first equation in (fTol l rewritten as (fl4] |. We now just 
need to check that ifTsl Assumptions 2.3 through 2.5] hold, so that we can apply ifTsl 
Theorem 2.4]. Notice that fl3, Assumption 2.3] is equivalent to Assumption 12.21 here 
and |13, Assumption 2.4] is satisfied due to Assumptions 12. ll here together with and 12.41 
Finally [13, Assumption 2.5] is satisfied thanks to our Assumption l2.4l Hence we get the 
desired rate of convergence. □ 

4. Approximations in cylindrical domains 

In this section we estimate the error arising in the localisation. To this end we prove a 
general result on the distribution of the exit time of diffusion processes form balls of some 
radius R. This, together with a comparison type theorem for the finite difference schemes 
from lfT3l allows us to estimate the localisation error 

Lemma 4.1. Let Wt be a dl -dimensional Wiener martingale. Let ^ be a measurable, 
valued random variable. Let (o't)tg[o,T] and (/3t)tG[o,T] be {^t)te[a.T] progressively 
measurable processes, where at is a d x d' dimensional matrix and j3t is a d' -dimensional 
vector such that \Ps\ds < oo andE, jaspds < oo. Let {xt)te[o,T] be given by 

dxt ^ I3tdt + (TtdWt, xo = ^■ 

Let Qs — ^(TsCr^. Let (•, •) denote the inner product in WK If for all t G [0, T] 

2xtl3t + IcTtp + {atxt,xt) < K{1 + jxtp) V - almost surely (15) 

and ifEe^ < oo, then there exists p. > 0, depending only on K and T, such that 

E sup e^^' < 3(l + Ee«'/2). (i6) 

te[o,T] 
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Proof. We begin by applying Ito's formula to the process xt- Thus 

dxl = 2xtatdWt + (2(a;t, A) + \at\^)dt 

and hence 

de"? =e"? [2xtcJtdWt + (2(a;t, pt) + Wt\'')dt + 2iatXt,Xt)dt] . 
Let ipt ■— exp{e^^*xf), where A > is a constant to be chosen later. Then 

dVt = e"^Vt [2xtatdWt + {2{xt,Pt) + P + {atXt,Xt))dt - Ax^dt] . 
Using (fTSl l we see that 

tpt<Tpo+ f 2e'^'ij,XsasdWs + [ e-^'iP^ {K{1 + x^) - \xl) ds. 
Jo Jq 

Let Tp denote the exist time of the process from a ball of radius R. Since stochastic 

integrals are local martingales we have 

.T' 

EIa^T' < ^lATpo + IEIa / e-^'ijs {K{1 + x^) - Xx^) ds =: El^Vo + / 

Jo 

for any stopping time T* < and any A E J^q. But letting i? — > oo we get the above 
inequality for any T* < T. We see that / is less then or equal to 

1; 



El, 



L{:.2<i}e~^"V's {Kil + a;2) - Ax^) ds + f l|^2>i}e-^^7/., {K{1 + x^J - Ax^) ds 

Jo 



Choose A > large such that, 2K < A. Then 

/ < EIa / l{a;2<i|e-^" exp(e-^"a;^) {K{1 + x^) - Aa;^) ds < El, 
Jo 

Hence for any stopping time T* < T we have 

ElAexp{e-^^'x^,) < E1a(1 + exp(C^)). 

Then, due to Gyongy and Krylov ifTTl Lemma 3.2], for any S E (0, 1) 

E sup exp(5e-^*a;2) < + Eexp((5C^)). 

te[o,T] 

Take S = 1/2. Let — exp(— Ar)/2. Then we see that 

E sup expifixf) < 3 ( 1 + Ecxp ( 
tG[o,r] V \2 

and this completes the proof. □ 

Corollary 4.2. LetT^ inf{i > : |a;t| > R}. If the assumptions of Lemma \4.1\ are 
satisfied then there exists fi > such that 

Proof. We will use Lemma l4n and Markov's inequality. 

P(r^ < T) = P I sup |xt| > i? I P I sup exp(^a;2) > exp(/ii?2) I 
\t€[0,T] J \telo,T] J 

<e-^^\ sup e^^* <3e-^^'(l + Ee«'/2)^ 
te[o,T] 

where /i > comes from l4.1l □ 
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Lemma 4.3. Let w be given by (|5]l, let g be a bounded Lipschitz continuous function of 
X and let Assumption \2J] be satisfied. Furthermore, let Ri > be given. Let qr^ be a 
function that is equal to g inside Bn-^, zero outside Bji-^^j^i and Lipschitz continuous. Let 

WR,{t,x):=¥.t,. sup (e-/r^K-")'i«g^J:cj,.)) , (17) 

T*^%[t.T\ ^ ' 

andletT*^^ inf{w G [i, T] : > Then 

\w{t, x) - WR, {t, x)\ < K^tAT'k < T}. 
Proof. Since the difference of supremums is less than the supremum of the difference, 

I:^\w{t,x)~ WR, (t, x) I < sup Et,, (ef"' "f"^-")'^" |g(a;T* ) - i^T- )\) ■ 

T'eilt.T] ^ ' 

Notice that for x E Br^, g{x) — qr^ [x) and that for all x E M"*, \g{x) — (a;)| < K. 
Hence 

/ < sup Et.,l{T- <T} f p(^.-^)duj^\ 
T'ei[t,T] 1 \ J 

Noting that p > concludes the proof. □ 

The following corollary is an immediate consequence of Lemma l43] and Corollary 14.21 

Corollary 4.4. Let w be given by ^ andwR-^ be given by (1171 ). Let R2 < Ri be a positive 
real number Then there exists /i > 0, such that for all {t, x) e [0, T] x Br^, 

\wit,x)-WR,{t,x)\<Ce-^"i+''-/^ 

In other words, by cutting the function g outside the ball of radius i?i we have in- 
troduced an error in the American put price that is decreasing exponentially as the ball 
increases in radius. Now we will show that, if the discrete problem is solved with a payoff 
function, that is zero outside a ball of radius Ri and the grid is restricted to a ball of radius 
R > Ri, then the difference between this solution and the solution of the discrete problem 
on the whole grid with the original payoff function decreases exponentially with R and i?i . 

Lemma 4.5. Let gi and g2 be functions of x, such that 51 < 52 < 00. Let ui,U2 be 
defined on Mt ond such that uie^'^'^' and U2e~'^'^' are bounded for some /i > 0. Let 
C > and assume that 

S'^ ui + LhUi + C, g ~ ui > max 5^ U2 + LhU2, g — U2 on Q (18) 

and Ui < U2 on A4t \ Q- If h < 1 then there is a constant f depending only on K,di, fi 
such that 

ui <U2 + TC on Mt, (19) 
for T e (0, f ). Ifui, U2 are bounded on Mt then il9\ holds for any t. 



Lemma 1431 is a special case of |13, Lemma 3.9], while following corollary is just a 
special case of |[13] Corollary 3.11]. 

Corollary 4.6. IfwT.h is the solution of ( fTOb then 

\wtm\ < C+ sup \g{x)\. 

Lemma 4.7. Assume that g(x) = 0/or > R. Then, for some 7 € (0, 1), 
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Proof. Let 7 e (0, 1) be some constant to be chosen later. Let ^(i) be defined recursively 
by: 

e(T) = l; m=l-^i{t + rT(t)) for t<T. 

Take an arbitrary unit I e M"^. Let 77(2;) ~ exp(7(x, I)) and C, = ^77. We would like to 
apply Lemma l43] Observe that by Taylor's theorem 

A,.,A v{y) = Di v{y) + /, °1 '^(2/ + s^fc)!/^ - 

+ sup |D^^77(2/ + s;fc)|. 

As |Zfc| < and 7 e (0, 1), we obtain that 

Bl7j = 7n'{l,lk)' <C^V and D^^ r? < ^77?. 

Hence 

sup |D^,?7(y + s;fc)| < (7777(2;) sup \exp{slk,lk)\ < Cjr]{y), 

sG{-h,h) se(-h,h) 

provided h < K. Thus Ah^^g^ rj < C^-q. Similarly (Jft^.^j. "H ^ djrj. Furthermore, since 
S'^i^v) = '''^^(7 ~ 1)C'7' we can see that, for sufficiently small 7 

^(J^C + UC-O < (T-i(7-l) + C7)C<0onHT. 

Let Q = {it,x) e Xt : (a:,/) < -R}. Since -i? > {x,l) = |a;|cos6' > 

Q C {{t,x) E A4t ■ \x\ > R}. Recall that g{x) — for > R. Hence for any constant 

C > 0, 

{5l CC + L„CC - rCC + rg) < on Q. 

For {t, x) e Mt \ Q, either t <T and {x, I) > -R oxt = T and (x, > -|a;|. Hence 
either 

C{t,x) > e''-'^oi-C{t,x) > e^'^l^l. 

In the first case, we know from Corollarv l4.6| that Wr.h is bounded by a constnat. Taking C 
large enough, e^^CC, > Wr^h- In the second case, t ~ T and so WT,h — 9- As g{x) — for 
|x| > i? we only need to consider < R and so for large C, e'^CC, > g = Wr,h- Either 
way, for C large enough, Ce'^C, > Wr,h on A4t \ Q- By Lemma l43] w^ >, < Ce^^C, in 
M.T- Since the choice of the unit vector I was arbitrary we can see that in A^y 

WT,h < Cre''''^ exp(— 7|a;|). 

Analogous application of Lemma |43] would yield that in M.t 

Wr,h > -CTe'^'^ exp{-j\x\) 

and hence complete the proof. □ 

Lemma 4.8. Let Assumption \2J] and \Z2\ be satisfied. Let R > Ri > and consider gn-^, 
a bounded Lipschitz continuous function that is zero outside B^^. Let w^]^ be a function 
satisfying ( llOl l on Mt and such that w^\^ = gji^ on A^t \ A^T- Let w^'l^^ be the solution 
to ( 1101 ) on Qji with w^'l^^ = gu^ on Mt \ Qr- Then for some constant 7 G (0, 1), 



ERROR ESTIMATES FOR FINITE DIFFERENCE APPROXIMATIONS OF AMERICAN PUT OPTION PRICE 1 1 

Proof. From Lemma EJ] we know that for \x\ > R, \w^l\ < Ce'^^^^~^l Since w^;^ 
satisfies (fTOb on Mt and w^'^^ satisfies ( fTOl i on Q^, we have, on Qj^, 

I S-T R,Rx , T R-Ri B,Ri\ A f xT Bi , V Ri Ri\ 

max [S^ w^ '^ + LhW^ f^ , g - w^ '^ j = ^ = max [d^ w^ l + hhw^ l,g - w^ lj 
< max [s^iw^X - Ce-^^^-^)) - Ce-^^^^-"'^), g - {w^X " . 

For (t, x) e Mt \ Qr either (t, x) e [0, T) x n A^t and so 
w?X-Ce^'"^-'''<0-9n.-wff^ 

or t = T and so 

<k - Ce-(«-«) = 9R. - Ce^(^-«) < = <f • 
Hence by Lemma l431 applied to ]^ — 

Pg7(fli-fl) andu;^f^ onQ^^ 
<1 < <f ^ + Ce^(^^-^) on >(t. 
Similar argument for + Ce'^^^'^~^^ gives 

w^X + Ce-^^^^-^) > wff' on A^t 
and so completes the proof. □ 

Proof of Theorem \2. 71 By Lemma [373l ( fTOl ) has a unique solution. We consider the optimal 
stopping problem wr^ given by (|5]), but with g replaced by ga^. By Corollary 14.41 there 
exists /i, > such that 

\w - WR, \ < Ce-''^i+^^/2 on [0,r] x Bb,- 

Let w^Jj be the solution of ( fTOb on A^t with g replaced by g^^j . This is the solution of the 
finite difference problem with g cut off outside a ball, but on a grid on the whole space. By 
Theorem l3.4l we know that 

\WR, - w^X\ ^ C-lri/* + /ii/2) on [q^t] x R". 

Finally, we consider wf'^'^, which is the solution of ( fTOl i on Qn and with g replaced by 
g/jj^. By Lemma l4!8] there exists 7 e (0, 1), such that 

l<k - </fM < on [0,T] x Br,. 

Hence, chaining the above three inequalities, we obtain 

\w - <f ^ I < C + + e^(«i-«) + e-Mfl?+fl|/2^ ^ 

which is the estimate for the error caused by both discretisation and localisation. □ 
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